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Abstract. According to standard models supernovae produce radioactive 44 Ti, which should be visible in gamma-rays follow- 
ing decay to 44 Ca for a few centuries. 44 Ti production is believed to be the source of cosmic 44 Ca, whose abundance is well 
established. Yet, gamma-ray telescopes have not seen the expected young remnants of core collapse events. The 44 Ti mean life 
of 89 y and the Galactic supernova rate of ^ 3/100 y imply ^ several detectable 44 Ti gamma-ray sources, but only one is 
clearly seen, the 340-year-old Cas A SNR. Furthermore, supernovae which produce much 44 Ti are expected to occur primarily 
in the inner part of the Galaxy, where young massive stars are most abundant. Because the Galaxy is transparent to gamma-rays, 
this should be the dominant location of expected gamma-ray sources. Yet the Cas A SNR as the only one source is located far 
from the inner Galaxy (at longitude 112°). We evaluate the surprising absence of detectable supernovae from the past three 
centuries. We discuss whether our understanding of SN explosions, their 44 Ti yields, their spatial distributions, and statistical 
arguments can be stretched so that this apparent disagreement may be accommodated within reasonable expectations, or if we 
have to revise some or all of the above aspects to bring expectations in agreement with the observations. We conclude that either 
core collapse supernovae have been improbably rare in the Galaxy during the past few centuries, or 44 Ti -producing supernovae 
are atypical supernovae. We also present a new argument based on 44 Ca/ 40 Ca ratios in mainstream SiC Stardust grains that may 
cast doubt on massive-He-cap Type I supernovae as the source of most galactic 44 Ca. 

Key words, nucleosynthesis; gamma rays: observations; supernovae: general; ISM: supernova remnants, abundances; 
astrochemistry-dust. 
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■ 1. Introduction scales. In this work we do this by asking if the current rate 
k> \ of supernovae in our Galaxy which produce radioactive 44 Ti 
r£ . Supernovae are the agents that drive the evolution of gaseous ifj ^ Hne whh expectations from other obs ervables and from 

■ regions of galaxies. As end points of the evolution of mas- supemova theory 

sive stars that have formed out of the interstellar gas, their 4 4 Ti decay offerfj a unique window to the study of supernova 

explosions eject matter enriched with freshly formed isotopes rates Spedfic aspects Qf M& window m . 
and stir interstellar gas. However, the explosions themselves 

are still not understood jBurrowsl feoOfi Ijanka et all I2003I) . 1- Gamma rays penetrate the entire galactic disk with little 

Parametric descriptions are used to describe the core collapses extinction 

("cc-SN", supernovae of types II and Ib/c) JWooslev & WeaveJ 2. 44 Ti gamma-rays reflect the current rate of supernovae, with 
ll995UThielemann et aiTl 996) as well as thermo nuclear explo- the 44 Ti mean decay time scale of r = 89 years; this present- 
sions of white dwarfs (supernovae of type la; No moto et aD day snapshot which has not yet fed back into chemical evo- 
Il997l) . A prominent issue in astrophysics is whether the super- lution can be directly related to the observable current pop- 
nova explosion itself is a well-regulated, robust physical pro- ulation of massive stars 

cess, or if intrinsic variability over a wider range of physical 3. Most 44 Ti is co-produced with 56 Ni in relatively frequent 

conditions are rather common. core collapse supernovae (cc-SN). The radioactive energy 

Supernova homogeneity by type can be studied in different of 56 Ni is responsible for well-observed supernova light 
ways. One approach is to analyze the rate of supernovae of a 4 - Nucleosynthesis of ^Ti is primarily from a-rich freeze-out 
specific type in different environments and over different time of nuclear statistical equilibrium and secondarily from sili- 
con burning 

Send offprint requests to: L.-S. The 5. The origin of abundant cosmic 44 Ca occurs mainly through 

Correspondence to: tlihsin@clemson.edu 44 Ti nucleosynthesis 
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6. Ti traces have been found in pre-solar grains which have 
been attributed to condensation within core-collapse super- 
novae 

In this paper we estimate what the gamma-ray sky of 44 Ti 
sources would be expected to look like by adopting an average 
AA Ti source model having a characteristic source event recur- 
rence rate, AA Ti yield per event, and spatial distribution. We 
compare this to the present-day gamma-ray survey and find 
apparent and serious conflicts. Then we analyze whether de- 
viations from these average expectations can occur from the 
known or expected variability of models and parameters in- 
volved. We use a Monte Carlo simulation of the expected sky 
image within reasonable distributions of parameters for that 
purpose. This leads us to discuss each of the relevant param- 
eters, which may explain an anomalous 44 Ti sky; these are, 
specifically: 

- the gamma-ray survey quality 

- statistical effects of small samples 

- the adopted supernova rates 

- the supernova explosion and nucleosynthesis models 

- the spatial distribution of supernova events 

- the deviations from smooth chemical evolution 

- the supernova origin of Galactic and solar ^Ca 

2. 44 Ti and Supernovae in the Galaxy 

Supernova rates in the Galaxy can be inferred from differ- 
ent observables. But observational incompleteness and bias re- 
quires that several assumptions are an essential part of such 
inferences. We compare here the observed 44 Ti sky with expec- 
tations for the occurrences of young 44 Ti emitting supernova 
remnants, as they result from astrophysical assessments of su- 
pernova characteristics for our Galaxy. We choose a time inter- 
val unit of 100 years for comparison of chemical history with 
supernova event rates. The extrapolation of the past history and 
yields of 44 Ti -ejecting supernova events should give a produc- 
tion rate which can be compared with recent supernova rates, 
as no other source of 44 Ca has yet been identified. 

Have the 44 Ca-producing supernovae been typical; is all 
of 44 Ca produced through radioactive 44 Ti; or do exceptional 
events contribute most of the 44 Ti ? Gamma-ray surveys for 
44 Ti sources, and presolar supernova grains provide ways to 
address this question. 

In the following Section we will then examine each of the 
critical assumptions in more detail. 

For a 44 Ti sky reference, we adopt the result from 
COMPTEL's survey in the 1.157 MeV band, which is the 
most complete survey to date JPupraz et all 1 9971 llvudin et alJ 
Il999i) . In this survey, one object has been clearly detected 
(340-year-old Cas A at a distance of 3.4 kpc), candidates 
at lower significance have been di scussed (most pro minently 
GRO J0852-4642 in the Vela region. llvudin et alJJ 19991)), and a 
weak signal from the Per OB2 association ( Dupraz et al. 119971) . 
Apparently, no bright young 44 Ti emitting supernova remnants 
are found in the inner region of the Galaxy (see Fig.Q. 

What do we expect the Galaxy to look like in 44 Ti emis- 
sion? 



2.1. 44 Ti from Supernovae 

Supernova observations directly demonstrate that these events 
are at the origin of 44 Ti production: The 1.157 MeV y-ray 
line following 44 Ti decay has been detect ed in the 340-year 
old Galactic s upernova remnant Cas A Jlvudin et all 1 1994 
IVinketalJ200ll) . Furthermore, SN1987As late light curve, ob- 
served in unique detail over more than 15 years, appears pow- 
ered by a similar amount of ^Ti (0.2-2.0 10~ 4 M o ), from mod- 
eling of radio active energy depo s ition and photon transport 
in the SNR dWooslev etal.lll989t iFransson & Kozmall2002h . 
Gamma-ray detection and proof of this interpretation is still 
lacking, but INTEGRAL'S recent observations may prove suf- 
ficiently sensitive. 

Presolar grains have been identified in meteoritic sam- 
ples through their unusual isotopic abundance patterns, and 
hold rich isotopic abundance d etail for c haracteri zing their con- 
densation environments (e.g. IClavton & Nittlerll2004l Izinnerl 
Il998h . SiC grains of the X-type are attributed to core-collapse 
supernovae from th eir large excesses in characteristic isotopes 
dNitfler et al.lll996l) . ^Ti -produced overabundance of ^Ca is 
found in these SiC X grains, indicating presence of 44 Ti at 
their time of condensation. Measured ^Ca/^Ca ratios are larg e 
dNittler etal.lll996t IClavton et al.lll997at iHoppe et alJhoOOh . 
This proves their supernova origin, on one hand, and likewise 
it proves that dust-producing supernovae may eject 44 Ti in sig- 
nificant amounts. Their attributio n to core collapse s upernovae 
rather than supernovae of type la (Nittler et al. 1996) is plausi- 
ble, if we believe that core collapse supernovae probably dom- 
inate 44 Ti production and because there is no direct evide nce o f 
dust condensation in thermonuclear supernovae jDweklF9 98i. 
It is likely that all of the observed X grains sample both differ- 
ent condensation environments and different production events. 
Therefore, their 44 Ca abundances cannot be interpreted in ab- 
solute terms as a measure of the mass of 44 Ti ejected per super- 
nova. 

The most plausible cosmic environment for production of 
^Ti is the a-rich freeze-out from high-tem perature burn- 
ing near Nuclear S tatistical Equilibrium ("e.g. IWooslev et all 
Il973t lArnettl ll996). The required high values for the entropy 
are found in core-collapse supernovae. Therefore the simplest 
plausible assumption is that core-collapse supernovae are re- 
sponsible for any substantial sources of 44 Ti. 

It is likewise plausible that traditional SNIa do not add 
significantly to nucleosynthesis at mass A=44, specifically no 
44 Ti, since their NSE fr eeze-out c onditions will not be favor- 
able for 44 Ti production (Timmes et alll995l) . We consider the 
symbiotic-star scenario as a rare subclass of SNI a, even tho ugh 
their 44 Ti yields may be large (iWooslev & Weaver! 1 19941) : so 
they indeed would be rare outliers, rather than typical 44 Ti pro- 
ducing supernovae. We will discuss uncertainties in each of 
these sites in more detail below (see !3.5> . 

2.2. Ca from Supernovae 

Measurements of cosmic isotopic abundances can be con- 
verted into isotope production rates; however, for all long- 
lived isotopes models of chemical evolution have to be applied 
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(e.g.lClavtorJl988tlPagell997tlTimmes et alll995l lMatteucci 
How can we use these for predicting the shortlived 44 Ti 
source appearance in gamma-rays? Its lifetime is too short for 
mean chemical evolution arguments. 

Over the time scale of chemical evolution of the Galaxy, the 
cumulative and averaging effects of different explosion types 
would integrate to a smooth pattern of "standard abundances", 
as they are observed throughout the universe. Evolutionary 
pathways for individual elements can be associated with the 
evolution of metallicity and the rates of different supernova 
types. The abundance of 44 Ca in solar-system matter thus can 
be translated through models of chemical evolution into a cur- 
rent average production rate of ^ 3 10 -4 M o of 44 Ti per 100 y 
ILeising & S hare 1990, see below). 

Chemical evolution calcu lations famines et al Jl 19951) us- 
ing computed ejecta masses IWoosley & Weaverlll995l) were 
shown to account reasonably for most solar abundances, in- 
cluding specifically 40 Ca, but failing by a factor of three for 
44 Ca; this fact is also evident from the ratio of 44 Ti to Fe in cc- 
SN models. But over short time scales where only a limited 
number of source events contribute, great variability among 
sources would present a significant difference of the present- 
day picture from the average. Therefore we may compare the 
expectations from long-term averaged 44 Ti production to the 
present-day 44 Ti source record imaged in gamma-rays, allow- 
ing for the short-term fluctuations with Monte-Carlo realiza- 
tions of the sources. 

The solar-system 44 Ca abundance is rooted p lausibly in 
Ti nucleosynthesis JWooslev et al 11 973Uclavtonl 1 982).' and 
is a result of the integrated Galactic nucleosynthesis prior to 
formation of the solar system. Given a time-dependence of the 
production rate, it can be normalized to the rate required to give 
precisely the measured solar abundance 4.5 Gyr ago, which 
also fixes the long-term average production rate today. Models 
of galactic chemical evolution, constrained by a number of ob- 
servables, provide us with that time dependence, albeit subject 
to a number of assumptions and parameter choices. Adjusting 
the 44 Ti yield in the aforementioned calculation JTimmes et alJ 
1 19951) upward by a factor of three to achieve the solar abun- 
dance, we infer a production rate of 44 Ca of about 3.6 10 -6 M 
yr -1 . Other considerations using different chemical evolution 
models (see Appendix O lead from the 44 Ca abundance to a 
current 44 Ti production rate: p( 44 Ti)=5.5 1O -6 M yr -1 , with 
full uncertainty range 1.2-12 1O -6 M yr -1 (see AppendixO. 

2.3. Supernova Rates 

Direct supernova rate measures have been made through cor- 
relations between supernova acti vity and other t r acers of the 
massive star content of a galaxy. Ivan den Berghl dl 99 ll) finds 
(2.62 +/- 0.8) h 2 m SN century" 1 . Forhi 00 = 0.75, the rate is 1.5 
+/- 0.8 SN century -1 . This rate is based on a combined study 
of galactic sup ernova remnants, historical SNe, and novae in 
M31 and M33. Cappellaro et al. ( 1993) refer to this rate as the 
best estimate. Ivan d en Bergh & Tam manJ dl 99 ll) find SNR = 
4.0 SN century -1 . The authors review supernova rates in ex- 
ternal galaxies and derive a specific supernova frequency, in 



units of 1 SNu = one SN per century per 10 ll) L (B), for var- 
ious galaxy types. If one assumes that the Galaxy is interme- 
diate between types Sab-Sb and types Sbc-Sd, the specific rate 
is 3 hj 00 SNu. For a Galactic blue-band luminosity of L(B) 
= 2.3xl0 10 L (B) (their Table 11) and h 100 = 0.75 we infer 
SNR = 4.0 SN century -1 . Their review paper also discusses 
estimates from internal tracers in the Milky Way: From radio 
supernova remnant (RSNR) statistics they infer SNR = 3.3 +/- 
2.0 SN century -1 . From the historic record of nearby (< a few 
kpc) supernovae in the past millennium they find SNR = 5.8 +/- 
2.4 SN century -1 . The large extinction corrections in the galac- 
tic plane make this small sample highly incomplete, which re- 
sults in large uncertainties in extrapolations to the full galac- 
tic disk. The authors also review efforts based on the pulsar 
birth rate, but extensive observational selection effects in com- 
bination with the strong and poorly understood evolution of 
lumino sity and beaming geometry (se e lLvn e & Graham-Smith 
Jl998l) and lLorimer & Kramerl (120041) 1 renders this method im- 
practical for estimating the galactic SNR. Continuing the stud- 
ies of van den Bergh and Tammann, Capp ellaro et alJ ill 9931) 
find SNR =1.4 +/- 0.9 SN century -1 , when scaling to external 
galaxies of similar type. The sample is obtained from surveys 
carried out at the Asiago and Sternberg Observatories. The au- 
thors provide an extensive discussion of the uncertainties of this 
method, which can exceed 200% for som e late type galaxies. 
More recently, Ivan den Bergh & McClurel JH94) find SNR = 
(2.4-2.7) h2. events per century. This estimate is based on re- 
evaluation of the extra-galactic SN rates obtained from Evans' 
1980-1988 observations. This method depends on a somewhat 
uncertain type of the Galaxy and the value of its blue-band lu- 
minosity, while the uncertainty due to the Hubble constant is 
now very small. Given the error analysis in the paper, the rate 
is uncertain by at least 30%. 

van den Bergh & McClurel ( 1 1994 in studying the super- 
nova rates of local spiral galaxies of types Sab-Sd of Evans' ob- 
servations estimated that 80%-90% of supernova i n that galax- 
ies are of Types Ibc and II. Recently Cappellaro (2003) com- 
bining five SN searches to include 137 SNe in 9,346 galax- 
ies estimates the SN type ratios in the Galaxy to be la : Ib/c 
: II = 0.22 : 0.11 : 0.67. From these observations, one in- 
fers a ratio of core-collapse to thermonuclear supernova of 
R=(II+Ibc)/Ia = 3.5. However, note that the Galactic histor- 
ical record in the last millennium shown in Table IB. II con- 
tains only two Type la SNR out of six SNRs. An often 
used alternative distrib ution over t ypes is la : Ib/c : II ' = 0.1 
: 0.15 : 0.75 (lHatano et alJll997l IPawson & Johnson! 1 19941 
ITammann et al.lll994t IHartmann et al.lll993l) . which implies a 
three times higher cc-SN fraction. In this work, we adopt this 
set of parameters, and note that the small 44 Ti yield of Type la 
renders our results insensitive to this ratio. 

Over the last millennium, the historic record contains six 
events (see Appendix IB.4I . which implies a rather low rate 
at face value. However Galactic extinction at visible wave- 
lengths and embedded supernovae will lead to large occulta- 
tion bias, and with extinction models plus Monte Carlo sim- 
ulations this historic record can be assessed to approximately 
agree with extragalactic rate determinations (see Section[3]and 
Appendix IA.2> . An often cited rate of galactic cc-SN of three 
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per century is consistent both w ith astronomical arguments 
(van den Bereh & McClure 1990) and with the rate inferred by 
Tim mes et alJ i 19951) from their chemical evolution model that 
produces solar abundances successfully. We adopt a supernova 
recurrence rate of 30 years as a baseline for our 44 Ti sky expec- 
tations. 
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Fig. 1. Maximum-entropy map of the Galactic plane (within latitude +30°) in the 44 7i' energy window ( 1 .066-1 .246 MeV) for the 
combination ofCGRO observ ations from 0.1 to 617.1. Two crosses mark the positions ofCas A and GROJ0852-4642. This figure 
is adopted from IvudTn] ( T/9991) . 
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Fig. 2. The expected Ti sky of 10 simulated galaxies of model A where the supernova recurrence time is taken to be 30 years and 
the supernovae ratio ofIa:Ib:II — 0.10:0.15:0.75. Simulating a 10 5 galaxy sky, a gamma-ray detector with a detection limit of 
\0~ 5 cm~ 2 s~ x would have a probability of detecting 0, 1, and 2 ^Ti sources of 0.0017, 0.012, and 0.037, respectively. A slightly 
better instrument than the 10 -5 cm -2 i _I detection limit would detect several 44 Ti sources. 
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2.4. Supernova Locations 

The 1.157 MeV y-ray line following 44 Ti decay has been de- 
tec ted_jn_lhe34£^y£arold_Gal^ remnant Cas 
A Jlvudin et alJll994t Ivhik et alJ|200ll) . COMPTEL's survey 
JPupraz et^llTl997rilvu^m^t^nfl99^) has resulted in other 
candidate sources, such as the so-called Vela junior SNR 
Jlvudin etalJ [l998). INTEGRAL'S inner-Galaxy survey has 
been studied with IBIS Ima ger data, which did not reveal 
a new source in this region (Renau d et alJl2004l) . The diffi- 
culties of MeV observations have thus not led to convinc- 
ing new 44 Ti rich supernova remnants, especially in the in- 
ner Galaxy region (l=0°+30°) where observations are deepest. 
Yet, Cas A seems an est ablished 44 Ti dete ction in C OMPTEL 
JSchonfelder et alJ2000l) and Be ppo-Sax JVinketal]200ll) and 
INTEGRAL/IBI S JVinkl 120051) measurements, while OSSE 
jThe et alJll99fj) and RXTE jRofhschild & LingenfelteJ2003h 
measurements of Cas A were not sufficiently sensitive. There 
exist no well-understood supernova remnants other than Cas A 
where the 44 Ti production issue can be tested. It is apparently 
the only supernova whose yield, age, and nearness makes 44 Ti 
visible in gammas. Naturally we ask ourselves if Cas A is a 
typical supernova or an anomalous case of a high- 44 Ti yield 
supernova? 

Current supernova models predict an amount of 44 Ti which 
is of the same order than what these observations suggest, 
though generally slightly less. Can we take this as a satisfactory 
confirmation of our understanding of core collapse supernova 
44 Ti production? These two identified core collapse events and 
their association with 44 Ti production appear to be in line with 
models which attribute 44 Ti production to the more frequent 
standard core collapse events but not to standard thermonuclear 
supernovae. Is this correct? 

If true, the location of 44 Ti sources should match the lo- 
cations of young massive stars which have rather short life- 
times. There is substantial evidence that massive star forma- 
tion occu rs in spiral arms an d predominantly in the inner 
Galax y JElmegreen et alJl2003UScoville et alJl200lUKennicutll 
1998). M assive stars can b e observed directly in the infrared 
(e.g. Maeder & Conti 1994), though extinction corrections are 
large in regions of dense clouds. Possibly an even better 
(though more indir ect) massive-star census can be derived from 
26 Al decay y-rays JDiehl et al][l9 95: Prantzos & Diehl 1996; 
lKnodlsedeJ2000l:lDiehl et alJ2005h . 26 Al is understood to orig- 
inate predominantly from massive stars and y-rays easily pen- 
etrate even dark clouds in star forming regions. So, do we see 
the 44 Ti sources in regions where we expect them to occur? Or 
do other factors which are not yet understood conspire to make 
44 Ti ejection a phenomenon of core collapse events occurring 
in special regions and environments? 

2.5. Estimating the 44 77 Sky Appearance 

If we want to estimate how the 44 Ti sky should appear in a 
gamma-ray survey, we need to follow a statistical approach, 
due to the rare occurrence of supernovae. We therefore apply a 
Monte Carlo approach of sampling plausible probability distri- 
butions for supernova rates, their 44 Ti yields, and their Galactic 



distribution, thus calculating a large statistical sample of possi- 
ble appearances of the 44 Ti sky. 

- The mean rate of supernovae in our Galaxy is taken as one 
event in 30 years. 

- 44 Ti producing supernova events can be of Type la, Ib/c, or 
II SNe. For this study, we choose a super nova type ratio, 
la : Ib/c : II = 0.10 : 0.15 : 0.75 JPawson & Johnson! 19941: 
lHatano et al.ll99l). 

- Supernovae are then assumed to produce 44 Ti according to 
models of that type. In our models, Type II supernova yields 
are generated according to stellar mass, which is drawn 
from a Salpeter-type initial mass distribution for M > 8 
M0; the ^Ti yiel ds per mass of the star are taken from 
Table 1 of lTimmes etal] dl99fil) . A typical yield is 3xl0~ 5 
M0 for M=25 Mq, but variations with mass are about a 
factor 2. The 44 Ti yields of Type lb SNe are uniformly dis- 
tributed between 3xlO -5 Mm and 9xl0~ 5 Mp with typi- 
cal values of 6xl0~ 5 M© dTimmes et al.l fl996). The 44 Ti 
yields of Type la SNe are uniformly distributed between 
8.7xl0~ 6 M and 2.7xl0~ 5 M so this range covers the 
44 Ti yield in the deflagration W7 model, the delayed deto- 
na tion WDD2 model, an d the late detonation W7DT model 
of BSTomoto et alJ (|l997). But instead of adopting the ^Ti 
ejecta masses directly from these models, we uniformly 
increase these by a factor 3 in order that chemical evolu- 
tion calculations reproduce the known solar 44 Ca/ 4() Ca ra- 
tio. This assumes that 44 Ti and 44 Ca nucleosynthesis are di- 
rectly related and that some unknown physics factor causes 
the computed 44 Ti masses to be uniformly low by a factor 
three. Note that 44 Ti yields and supernova rate determine 
the brightness scale of our 44 Ti sky expectations. 

- For each supernova type, we then adopt a parent spatial dis- 
tribution, from which we draw random samples to deter- 
mine the location of the event. The spatial distribution for 
core-collapse supernovae (Types Ib/c and II) is taken, from 
either an exponential disk with scale radii between 3.5 and 
5 kpc, or from a Gaussian-shaped ring at 3.7 kpc radius 
with a thickness of 1.3 kpc, both placing young stars in the 
inner region of the Galaxy as suggested by many observ- 
ables. Thermonuclear supernovae (Type la) are assumed to 
arise from an old stellar population, hence we adopt the 
distribution of novae as our parent Galactic distribution for 
those rare events . We adopt the composit e disk-spheroid 
nova model from Higdon & Fowler] 41987) (See Appendix 
|A]for more detail on the variety of modeled distributions). 

- We then randomly choose the individual supernova age 
within the last millennium, for deriving its 44 Ti decay 
gamma-ray brightness. This assumes that the factors gov- 
erning the mean recurrence rate have not changed in the 
past 10 3 yr. 

Randomly selected sample results from 10 5 such realiza- 
tions of an expected 44 Ti sky are shown in Fig. |2] to be com- 
pared with the observed 1.157 MeV COMPTEL image shown 
in Fig. [jlvudin et al. 1999). We illustrate how typical these ex- 
pected images would be, by showing the distribution of source 
brightnesses over a much larger number of Monte Carlo sam- 
ples (Fig. |3] left): A ^Ti flux above a representative limit of 
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Fig. 3. Cumulative ^Ti gamma-ray line flux distribution of supernovae with f y >lxl0 12 cmT 2 s~ l according to our Model A (see 
Appendix \A.2\ with supernovae type ratio ofla.Ib.II — 0.1 :0. 15:0.75 (left). Expected average number of supernovae per sample 
galaxy with their 44 Ti gamma-line fluxes above detection limits, for a supernova recurrence time of 30 yrs, irrespective of their 
position in the sky ( right). The detection limit of COMPTEL instrument is shown as a dotted line. 



10~ 5 ph cm _2 s _1 occurs in 19% of our 44 Ti gamma-line flux 
distribution. The expected number of 44 Ti point sources per 
galaxy lies above the gamma-ray survey sensitivity limit, i.e., 
we do expect typically 5-6 positive detections of 44 Ti sources 
(see Fig. fright). 
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Fig. 4. The longitude distribution of supernovae for four Ti 
gamma-ray line flux bandwidths in Model A of Appendix 
\A.2\ with supernovae type ratio of la: lb: II — 0.1 :0. 15:0.75. 
The equal-bin width is 20 degrees in longitude. Only the 
f y <lxl0~ 6 cmT 2 s~ l distribution is normalized and multiplied 
by a factor of 0.298. 

It is immediately evident from Figures|2]and[3]that expecta- 
tions based on such seemingly plausible assumptions look very 
different than the observed 44 Ti sky: Fig.|2]shows ~ 4-7 observ- 
able 44 Ti sources in an area of the Galaxy that contained none 
(see Fig.0 above the observable flux limit used for compari- 
son. The brightest source in that realization (large-filled circles) 
has l y ~ 10 _4 cm _2 s _1 , which would have been seen as a ~10 cr 
source by COMPTEL, and would have been detected already in 



INTEGRAL'S inner-Galaxy survey JVinkll2004l> . The majority 
of our calculated samples lead to this same type of conflict (see 
Fig. |3- The probability of having no sources within the central 
galaxy is small for current surveys, as illustrated in Fig.[4](for 
model A): For a survey down to \ y ~ 10~ 5 cirT 2 s~', ~12% of su- 
pernovae appear within longitudes 0+60°, whereas only 72% of 
supernovae is within that volume . An interesting feature, how- 
ever, is that the longitude distribution for f y <lxl0~ 5 cm" 2 s _1 
is dominated by sources from Type la (the bulge) and farther 
away SNe, while the distribution for f y >lxl0~ 5 cm _2 s _1 is 
dominated by sources from Type lb and Type II (disk) and 
nearby SNe. This feature also can be seen in Fig. [3] where 
around f y ^lxlCT 5 cm~ 2 s -1 the distribution changes its slope. 

It is clear that something is wrong with at least one of these 
assumptions: 

- the rate of core collapse supernovae in the Galaxy is 3/100 y 

- a core collapse supernova produces ^1O _4 M of 44 Ti 

- supernovae of each given type sample a relatively narrow 
range of the 44 Ti yield, i.e., we do not consider extreme 
outliers event (nor do we consider a rare but superhigh yield 
source class). 

- the locations of core collapse supernovae is axisymmetric 
i.e., we assume star formation in spiral arms is not domi- 
nant. 

In the following we examine these questions, and seek pos- 
sible explanations for the apparent conflict. 

3. Resolving the 44 Ti Sky Conflict 

3.1. Gamma-ray Observations 

The 44 Ti gamma-ray sky can be studied in lines from the pri- 
mary decay to 44 Sc at 67.9 and 78.4 keV and from the following 
decay to 44 Ca at 1.157 MeV. In this latter line, the COMPTEL 
imaging tele scope had reported the pioneering detection of 44 Ti 
from Cas A dlvudin et all 19941) . clearly showing a point source 
image at 1.157 MeV as well as the line in a spectrum from 
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this source. This detection had initially created some contro- 
versy b ecause other gamma-ray in struments apparently did not 
see it dThe et alJll99fl iRothschild et alJ ll998). We now be- 
lieve that this is due to the high initially-reported COMPTEL 
gamma-ray flux value, reduc ed later with better statistical ac- 
curac y of the measurement (Iv udin et al.lll999l iDupraz et al 



2001 



1997). The independent detections with BeppoSax ( Vink et al 



and with INTEGRAL/IBIS dVinkll2005l) in both lower- 
energy lines (Beppo-SAX) and in the 67.9 keV line (IBIS) 
now consolidate the 44 Ti detection from Cas A, but also sug- 
gest that indeed the 44 Ti gamma-ray flux of Cas A is in a 
range between 0.8 and 3.5 10 s ph cm _2 s _1 ; an "avera ge" o f 
(2.6+0.4±0.5) 10~ 5 ph cnrV 1 has been derived dVinkl2005l) . 
An upper limit from INTEGRAL/SPI reported from first stud- 
ies is consistent with this flux value, and may su ggest that the 
^Ti line is broader than 1000 km s" 1 JVinkl2005b . 

COMPTEL' s sky survey allowed for mapping of the plane 
of the Galaxy in the 44 Ti line (Dupraz et al. 1996, Iyudin et al. 
1999). Secondary features in these COMPTEL 44 Ti maps kept 
the discussion about statistical significances a nd systematic un- 
certainties alive ("see lSchonfelder et alJJ2000h for a comparison 
of Cas A to RXJ0852, a promising second source candidate, 
llvudin et all Jl998hL The COMPTEL point source detection 
algorithm dde Boer et ail 1 19921) has been tested with simula- 
tions over the full sky: Likelihood statistics has been verified 
to reproduce the expected number of artificial sources for a full 
sky survey, as the noise level is approached. The problem is 
that in the range of all the 44 Ti gamma-ray lines, all gamma- 
ray telescopes suffer from a large background from local ra- 
dioactivity hiduc£d_j2y_«)smi£j^ 

strum ents dGehrels & Michelsonll 19991 IWeidenspointner et alJ 
120021) . Determination of this background is crucial. For an 
imaging instrument, this can be done rather well by interpo- 
lation of imaging signatures from adjacent energies. But fur- 
thermore, for COMPTEL the 1.157 MeV line of 44 Ti is not 
far above its lower energy threshold, and in fact imaging se- 
lections strongly affect the sensitivity of the instrument up to 
1.5 MeV. Nevertheless, imaging analysis in adjacent energy 
bands should all experience similar problems, and therefore 
differences between images in the 44 Ti band and in neighbor- 
ing energy bands can be attributed to ^Ti rather than contin- 
uum sources or instrumental artifacts, once they are confirmed 
to be point-like sources (instrumental background lines would 
in general spread over data space and hence lead to extended or 
large-scale artifacts in the image, hence impact on the flux mea- 
surement rather than on point source detection). Detailed com- 
parisons of results for different energy bands, data subsets, and 
selections hav e led to the more cauti ous report about the Vela- 
region source ( Schonfel der et al.l20001) . so that still no convinc- 
ing second 44 Ti source clearly above the noise level is claimed. 
The INTEGRAL limits for 44 Ti lines from the Vela region 
source are now close t o the reported C OMPTEL flux value for 
this candidate source dRenaudll2004t Ivon Kienlin et alJEooH) 
and we therefore consider only one single 44 Ti source (Cas A) 
as being detected down to flux levels of 10~ 5 ph cm _2 s _I . 

For the low-energy lines measured by other instruments, in 
addition t o the instrumen tal background the underlying contin- 
uum (see The et al. ( 1999) for Cas A) presents a major source 



of uncertain ty (see lVink etaD d200lb : fvink & Laming! d2003l) : 
IVinkl 12005)1. This becomes even more of a problem for recom- 
bination line feature s from ^Sc at 4.1 keV, which have bee n 
studied with ASCA dTsunemi et al.ll2000t llvudin et al.ll2005h . 
yet without clear detections (though tantalizing hints have been 
discussed for the Vela-region source). 

For our quantitative comparisons of the 44 Ti supernova rate 
with the historical record of supernovae in the last millennium 
(see Appendix IB. 2> . we make use of the COMPTEL survey. 
We avoid the regions where reduced exposure might eventu- 
ally lead to increased artifact levels, and concentrate on the in- 
ner Galaxy where exposure for the COMPTEL sky survey is 
deep and homogeneous. We also avoid using data after the sec- 
ond re-boost of the satellite, whereaft er the activation level o f 
22 Na had increased substantially ( Weidenspointn eTet alJ20 02). 
This leaves us with a dataset covering the sky range of \l\ <90°, 
\b\ <30° using 7 years of the 9-year sky survey. In these 
gamma-ray data, the 44 Ti source sensitivity should be rather 
well-behaved and useful for our study. Note that the first inner- 
Galaxy s urvey from IN TEGRAL (longitudes of +20° around 
the GC) dRenaudl l2004) is consistent with our dataset in that 
also there no s ource is found to flux levels of 10 ph cm _2 s _1 
dRenaudl2005l) . 

3.2. Galactic Supernova Rate 

Many determinations of the galactic star formation rate (SFR) 
and supernova rate (SNR) have been made in past decades 
(see Diehl et al. 200 5 for a compilation of estimates, and 
Stahler & Palla (2005) on general astrophysical aspects of star 
formation). In this work we are concerned with the 44 Ti from 
supernovae, so we need the formation rate of massive stars 
(above, say, 10 M ). The conversion between SNR and SFR is 
a sensitive function of the Initial Mass Function (IMF), and val- 
ues given in the literature thus vary depending on the author's 
chosen IMF. For a generic transform ation e quation, we use the 
calibration from lMcKee & Williams! dl997l) : SFR = 196 SNR, 
where the SFR is measured in M yr -1 and the SNR in events 
per year. A star formation rate of 4 M Q yr -1 thus corresponds to 
a supernova rate of two events per century. These rates only in- 
clude core collapse supernovae (Type II and Ibc), but not SNIa. 
Generally these supernova rates are averages over time scales 
much longer than 44 Ti decay. 

Many papers discuss the star formation (rate) history (SFH) 
in relative terms, or the star formation rate surface density (M 
yr kpc~ 2 ) in the solar neighborhood and its radial depen- 
dence. None of these papers is useful (for our purpose) without 
an absolute calibration, based on a model of the galactic distri- 
bution of star formation. 

Generically, the SFR is obtained from a tracer that can be 
corrected for observational selection effects and is understood 
well enough so that possible evolutionary effects can be taken 
into account. One either deals with a class of residual objects, 
such as pulsars or supernova remnants, or with reprocessed 
light, such as free-free, H-alpha, or IR emission that follows 
from the ionization and heating of interstellar gas and its dust 
content in the vicinity of the hot and luminous stars. One must 
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be careful to include time-dependent effects, as the afterglow 
of an instantaneous starburst behaves differently than the steady 
state output from a region with continuous star formation. Here, 
we are concerned with an average star formation rate. 

To set the stage, let us collect SFR values from the litera- 
ture. ISmifh et alJ i ll 9781) concluded SFR = 5.3 M yr -1 , while 
TalbotNl98dTrinds a very low value SFR = 0.8 M G yr _1 , and 
Guesten & MezgeJ lll982h report a very high value of SFR = 
13.0 Mp yr~' Later papers appear to converge on a median 
value: iTurnel ill 984 finds SFR = 3.0 M yr" 1 , and iMezg eJ 
Jl987l) finds 5.1 M Q yr -1 . Measurements in the past decade 
confirm this moderate rate: McK3 dl989h derives SFR = 3.6 
M yr -1 from the analysis of thermal radio emission (free-free) 
from HII regions around massive stars. This emission is di- 
rectly proportional to the production rate of ionizing photons, 
which in turn is directly proportional to the SFR. He pointed 
out that this method is very sensitive to the slope of the high- 
mass IMF. We also note that the method depends on stellar 
model atmospheres in conjunction with models for massive 
stars, which change with treatments of mass loss, rotation, and 
convection. This paper also briefly discusses the use of the far- 
IR luminosity, due to warm dust heated by the absorption of 
photons from massive stars. The author uses the measured IR 
luminosity of the Galaxy of 4.7 xlO 9 L G (from iMezgeill 19871) 
to derive SFR = 2.4 M Q yr -1 . 

We must convert between star forma tion rate (SFR) and 
supernova rate (SNR) (see abovel. lMcKee & Williams! i 19971) 
prom ote the value SFR = 4.0 M yr -1 , and based on the Scalo's 
IMF JScalol ll986) convert this rate into a total number rate 
of 7.9 stars per year. They assume that all stars above 8 M 
become supernovae, corresponding to a supernova fraction of 
2.6x 10~ 3 . The mean stellar mass is < m > = 0.51 M . The 
corresponding cc-SN rate is 2 per century. This value is also 
supported by a completely independent method, based on the 
production of radioactive 26 Al in cc-SNe, whic h can be traced 
though its gamma-ray line at 1.809 MeV. Tim mes et al] 0997) 
use this method and conclude SFR = (5 +/- 4) M yr -1 , uti- 
lizing the 26 Al line flux derived from COMPTEL. The steady- 
state mass of 26 Al obtained in their work is in the range 0.7 - 2.8 
M , consistent with the value presented in the recent study by 
Diehl et al. (2005). Based on the S alpeter IMF in the range 0.1 
- 40 M and the 26 Al yields from lWooslev & Weaver! 0995) 
[which do not include contributions from the Wolf-Rayet wind 
phase] the authors derive the above quoted SFR and an asso- 
ciated cc-supernova rate of 3.4 +/- 2.8 per century. They ne- 
glect hydrostatically produced 26 Al that is injected into the ISM 
in massive star winds, which causes their SFR to be overesti- 
mated. The large uncertainty is mostly due to the steady-state 
mass of 26 Al inferred from the COMPTEL flux. INTEGRAL 
data presented by Diehl et al. (2005) have significantly reduced 
the error in this key quantity, and with wind yields included 
the latter study finds SNR = 1.9 + 1.1 supernovae per cen- 
tury, which corresponds to S FR - 3.8 + 2.2 M yr -1 , similar to 



the value given in McKee & Williams|ill997l). This SFR is very 
similar to the one obtained for M51 ( Calzetti et al.l l2005). and 
thus places the Galaxy in the group of quiescently star-forming 
galaxies. 



Th e most recent paper addressing this issue is by Reec 
(2005), who does not derive the SFR, but states that the galac- 
tic supernova rate is estimated as probably not less than 1 nor 
more than 2 per century. Using the conversion factors from 
McKee & Williams (1997), one infers that the SFR is in the 
range 2-4 M yr -1 . Reed uses a sample of a little over 400 03- 
B2 dwarfs within 1.5 kpc of the Sun to determine the birthrate 
of stars more massive than 10 M . The galaxy wide rate is 
derived from this local measurement by extrapolation based 
on models for the spatial distribution of stars, a model for 
galactic extinction (to accomplish corrections for stellar mag- 
nitudes), and a model of stellar life times. Reed emphasizes 
various sources of errors, such as lacking spectral classifica- 
tions of some bright OB stars, the (poorly known) inhomoge- 
neous spatial structure of extinction as well as stellar density, 
and non-unique connection between mass and spectral type. 
Finally, Reed also draws attention to the fact that one would 
have to include B3 dwarfs as well, if t he lower mass limi t for 
supernovae is 8 M and not 10 M (see He ger et al.l ( 120031) for 
comments on this mass limit). The OB-star catalog of the au- 
thor was used to perform a modified V/Vmax test to obtain a 
present-day star count as a function of absolute V-band magni- 
tude. From the stellar life times, and the assumption of steady 
state, the local birthrate follows. A double exponential model 
(in galactocentric radius and scale height above the plane) of 
the spatial distribution of these stars (which includes an inner 
hole of radius R = 4.25 kpc) ultimately leads to a total birthrate 
of 1.14 OB stars per century. Variations in the size of the hole 
change this number significantly, which leads the author to fi- 
nally claim a rate of 1-2 supernovae per century. 

In the context of our interpretation of the 44 Ti observa- 
tions, we would argue that Supernova rates between one and 
three cc-SNe per century are consistent with the large set of 
studies reviewed above. To solve the 44 Ti Sky conflict with 
a choice of the SNR (or SFR), an extremely low rate outside 
this range would of course explain the absence of 44 Ti gamma- 
ray line sources in the sky. However, chemical evolution argu- 
ments for 44 Ca would then require correspondingly higher ^Ti 
yields which are not supported by explosive nucleosynthesis 
studies (as discussed below) and which in any case would lead 
to higher fluxes from supernova remnants and thus again to a 
source count that exceeds the observed count of ~ 1 . The ob- 
servational constraints are on the product of rate and yield. The 
natural solution to the problem may be the very rare events with 
extremely high yields, which is discussed below. In that case it 
is of course totally unclear what to use for the spatial distribu- 
tion of these events, and "unusual" positions of a gamma-ray 
line source on the sky (such as that of Cas A) would be hard to 
interpret. 

But perhaps we have overlooked another option. We know 
that the galactic s tar formation process i s strongly correlated 
in space and time (Elm egreen et alJl2003l) . Could it be that the 
Galaxy just had a brief hiatus in its SFR? It would not take too 
much of a pause (say a few centuries) to explain the absence of 
bright 44 Ti sources if the past few centuries were very untypical 
with respect to the SFR (or SNR). This possibility is included 
in our Monte Carlo simulations, which show that this is not a 
likely solution when one simply considers Poisson fluctuations. 



10 



The et al.: 44 Ti Supernova Remnants 



This solution is thus not acceptable, unless one can point out a 
physical cause of the hiatus in the recent SFR. 

3.3. Galactic Supernovae Locations 

Likewise, we may wonder about the possibility of very large 
spatial fluctuations. Could it be that the recent Galaxy exhibits 
an average star formation, rate-wise, but a lopsided distribu- 
tion in space. If the opposite side of the Galaxy currently forms 
stars, and regions in the solar sector are relatively inactive (/ = 
90 - 270°), we expect to detect fewer gamma-ray line sources 
(because of their somewhat larger average distances). But at 
the same time the supernovae would also suffer from enhanced 
extinction, and matching the historic SN record would require 
a much higher rate. We have simulated the effects of a lopsided 
profile with a von Mises Distribu tion, the analog of a Ga ussian 
distribution for circular data (Mardia 1975; Fisher 1996). This 
function allows us to change from an axis-symmetric galactic 
distribution to a one-parameter distribution (measured by a pa- 
rameter k, where k=0 corresponds to the uniform circular dis- 
tribution) in a chosen direction which we choose it to be the 
longitude I = 0° direction. We re-simulated the 10 5 Galaxies 
sample as described in Sec. 12.51 and find that the detection 
probabilities (stated in the caption of Fig. [2} change to 0.0042, 
0.023, and 0.062 for parameter k=0.5 and 0.0068, 0.035, and 
0.087 for parameter k= 1 .0, respectively. The probability curves 
for 44 Ti y-line source detection (see Fig. IB.4I > are shifted to- 
wards higher rates, as expected, but the overall likelihood of 
these models decreases. The constraint from the historic record 
demands new rates that are even larger, as the extinction correc- 
tion affects the results more strongly than the D 2 distance ef- 
fect for the flux. The combination of these two constraints make 
lopsided models less acceptable than axis-symmetric ones. 

A lopsided model would also make Cas A even more spe- 
cial, regarding its unexpected location on the sky. To alleviate 
this problem we simulated lopsided star-forming galaxies in 
which the solar sector was the more active. The detection prob- 
abilities of 0, 1, and 2 44 Ti sources of this model are 0.0005, 
0.003, and 0.012 forparameterk=1.0, respectively. This shows 
that the number of the most probable 44 Ti y-line source detec- 
tion in the Galaxy is >2 detections and the model is less prob- 
able than the model used for Fig. |3 for consistency with the 
observed 44 Ti sky. 

Is it reasonable at all to consider one-sided star form- 
ing galaxies? That major merger events should be able to 
tidally induce a lopsided starburst activity is perhaps obvi- 
ous, but the Galaxy is not undergoing such an event. However, 
lRiidjrick^etaljj |2000l) have shown that even minor mergers 
Jlbata et all 19951) generically termed "weak interactions" may 
lead to a boost in the star formation rate correlated with their 
lopsidedness. Another mechanism for the creation of non- 
symmetric star formation patterns is the interaction between 
odd numbers of spiral d ensity waves, as it may be at work in 
M51 (He nry et all20 03). We do not advocate such an asymme- 
try for our Galaxy, but just wanted to consider this real possi- 
bility as one of the potential fixes for the 44 Ti sky problem. Our 
simulations indicate that even such an extreme solution does 



not work, as the various combined constraints operate against 
each other. While a lopsided Galaxy helps on the gamma-ray 
source count side, the historic record is harder to explain if re- 
cent supernovae are located preferentially on the far-side. 

3.4. Galactic Extinction Map 

Analyzing the consistency of the supernova rates derived from 
the Galactic historical record and the COMPTEL's gamma- 
ray map (Appendix IB.4I >. the consistency would improve if 
the SN rate from historical record is smaller. This could be 
realized if the true Galactic extinction map is lower than the 
visual extinction map we used here. Recently there are two 
published Galactic extinction maps that are useful for the 
type of study of this paper. The optical reddening model of 
Mendez & van Altenal (17998). which is based on Galactic dust 
distribution model, makes use of t he sa me optical sky sur- 
veys implemented by Hakki kTet alJ dl997l) in addition to some 
other restricted surveys. However this extinction map is reli- 
able only for solar neighborhood within 6 kpc. Another large 
scale three-dimensional model of Gal actic exti nction based on 
the Galactic dust distribution of Drim mel et alJ J2003I) has been 
shown to give a good agreement with the empirical extinction 
derived from NIR color-magnitude diagrams within 0.05 mag- 
nitude and furthermore it is reliable for a distance up to ^8 kpc. 
This extinc tion model gives a l arger magnitude of extinction 
than that o fHak kilaetal](ll997l) for longitude |Z| < 1.5°and for 
most pointing directions from the Sun for distance larger than 
6 kpc. For distance less tha n 5 kpc, the ex tinction of this map 
is smaller than that of Hakk ila et alJ {7997) which could give a 
better agreement between the SN rates of the historical record 
and of the COMPTEL's gamma-ray map. 

3.5. The SN model 
3.5.1. Lifetime of 44 Ti 

Although the lifetime of 44 Ti measured in laboratories had 
exhibited a large uncertainty since its first measurement in 
1965, a compilation of five recent experiments performed af- 
ter 1998 gives an averaged lifetime of 87 + 1 yr, where the 
quoted error is o f statistical and of one standard d eviation (see, 
e.g., F igure 5 of Hashi moto et alJ ( 120011) and also lGorres et alJ 
( 1998)). Apparently, this small uncertainty in the measured 44 Ti 
lifetime does not affect the discrepancy discussed here. 

It is noted that the above-mentioned lifetime measured in 
laboratories is for neutral atoms. Since 44 Ti is a pure orbital- 
electron-capture decay isotope, its lifetime depends on the elec- 
tronic environment in the evolutional course of a supernova 
remnant. For example, a fully-ionized 44 Ti is stable, and the 
lifetime of 44 Ti in the Hydrogen-like ionization state becomes 
long er by a factor of 2.25 than that of the neutral 44 Ti (see 
Moti zuki & Kumag ail (12004)1. Let us briefly consider the effect 
of 44 Ti ionization on our problem. 

In young supernova remnants, the reverse shock propagates 
inward through the ejecta and the resulting increase in temper- 
ature and density may lead to highly ionized ejecta material 
through thermal collisions with free electrons. A high-degree 
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of ionization may then result in a longer lifetime of Ti, which 
would significantly alter the inferred 44 Ti mass. In fact, H-like 
and H e-like Fe ions have been observed in Cas A (see, e.g., 
iHwang etalJ E004). Because the electron binding energies of 
Ti are smaller than those of Fe, it is easier to ionize Ti than 
Fe. Accordingly, 44 Ti atoms in Cas A may be expected to be in 
such high ionization states at least in part if they are accompa- 
nied by the highly ionized Fe (this is expected because 44 Ti is 
synthesized at the same location as where 56 Ni is also produced 
in the innermost region of a supernova). 

Since the present-day radioactivity was entirely affected by 
the history of various ionization stages and their duration time 
for which the 44 Ti has experienced through the evolution, de- 
taile d discussion requires n ume rical simula t ions a s was done 
by iMochizuki et all i 19991) and iMochizukU feOOll) . However, 
we can get a rough idea of the ionization effect on the radioac- 
tivit y by using the result of si mple linear analysis, i.e., eq. (7) 
of lMotizuki & Kumaeail J2004 : 

AFy/Fy = (t/r- 1)At/t, (1) 

where we have replaced the radi oactivity A and the d ecay rate 
A appeared in eq. (7) of lMotizuki & Kumagail (120041) with the 
y-ray flux F y and the 44 Ti lifetime t, respectively. In eq.([Q, t 
is the age of a SNR, AF y is the change of the flux by ionization, 
At is that of the lifetime. 

Note that At is always positive because the ionization al- 
ways increases its lifetime. As was pointed out by the above 
authors, the sign of AF y is then determined by that of the term 
in the parenthesis in the right-hand side of eq.(^. This means 
that the flux is enhanced by the ionization when a SNR is older 
than the 44 Ti lifetime, and that the flux is reduced when it is 
younger. 

Our concern here is whether the effect of ionization on the 
lifetime of 44 Ti in SNRs can reduce the disagreement between 
the observed 44 Ti Galactic map and the model's map or not. 
From the above arguments, we can easily understand that the 
discrepancy may be diminished if the y-line fluxes in Fig. |2] 
could be smaller which may be realized if most of the y-ray 
detected SNRs in Fig.[2]are younger than the 44 Ti lifetime. 

To get a rough idea, we performed a calculation in which all 
parameters are the same as employed for Fig.|2]except 1) the 
fluxes are multiplied by a factor 0.5 for SNRs with ages less 
than 100 y, and 2) the fluxes are multiplied by a factor 2 for 
SNRs with ages between 200 and 400 y. The selection of 200- 
400 years old SNRs as enhanced targets here is because the 
effect of the ionization due to the reverse shock is considered 
to be distinguished for these ages and the further inclusion of 
the enhanced-flux effect on SNRs olde r than 400 y only makes 
the discrepancy larger (see Mochizuki et al. ( 1999) for details). 
Simulating a 10 5 galaxy sky, we found that a y-ray detector 
with a detection limit of 1 x 10~ 5 ph cirT 2 s _1 would have a 
probability of detecting 0, 1, 2 44 Ti sources of 0.0012, 0.008, 
and 0.026, respectively. Therefore, from this simple analysis it 
is suggested that the disagreement cannot be compensated by 
the ionization effect; in effect it becomes worse in our simple 
calculations above than the reference calculations of Fig. [2] 

A more precise estimate requires the knowledge of the tem- 
perature and the density evolution of a supernova remnant, and 



the distribution of Ti in it. However, any detailed calculations 
taking into account the retardation of 44 Ti decay due to ion- 
ization will not alter the situation better: In any case, the older 
SNRs whose fluxes may be enhanced always dominate in num- 
ber the younger SNRs whose fluxes may be decreased. 

3.5.2. Nucleosynthesis reaction rates 

Estimates of yields of 44 Ti from nucleosynthesis in super- 
novae depend crucially on key nuclear reaction rates, and un- 
certainties in these rate s lim it our ability to constrain the su- 
pernova rate. iThe et all Jl998t) studied the sensitivity of 44 Ti 
yields in alpha-rich freezeouts to uncertainties in nuclear reac- 
tion rates. They did this by computing the alpha-rich freeze- 
out with reference values for the reaction rates and then com- 
paring these results with ones from calculations with individ- 
ual rates varied upwards and downwards by a factor of 100 
from their reference values. The results were that the produc- 
tion of 44 Ti was most sensitive to the rates for the follow- 
ing reactions: 44 Ti(a, pf 1 V, a(2a,y) n C, u Ti(a, y) 48 Cr, and 
45 V(p, y) 46 Cr for matter with equal numbers of neutrons and 
protons (77 = 0). For neutron excess 77 greater than zero, the 
importance of the reaction 45 V(p, y) 46 Cr drops, but other reac- 
tions become more important. In particular, these reactions are 
n C(a, y) ie O, 40 Ca(o<, y) 44 Ti, 27 Al(o:, «) 30 P, and 30 Si(a, «) 33 S. 

For our purposes, the relevant question is how much the 
44 Ti may vary from current supernova models given thes e 
uncertain ties. Motivat ed by the work of IThe et alJ dl998h . 
ISonzogni et alJ fcOOOt) measured the cross section for the 
44 Ti(a, pf 1 V reaction at the astrophysically relevant energies. 
They found that the experimental cross section for this reac- 
ti on was a factor of two l arger than i n the rate compilation 
of iThielemann et alJ Jl987h used in thelThe et alJ (fl998) cal- 
culations. From this result, ISonzogni et alJ J2000I) inferred a 
25% reduction in the amount of 44 Ti produced in alpha-rich 
freezeouts in supernovae. Other of the key reactions found by 
IThe et all lll998h had similar sensitivities of 44 Ti yield to re- 
action rates; therefore, if other experimental reaction rates are 
also a factor of ~ 2 different from the theoretical calculations, 
we can expect similar ~ 25% effects on the 44 Ti yield. From 
these results, we might thus conservatively expect the 44 Ti yield 
to be uncertain by less than a factor of ~ 2 due to reaction rate 
uncertainties. Such a conclusion is supported by the study of 
the reaction rate sensiti vity of nucleosynthesi s yields in core- 
collapse supernovae by Hoffman et alJ dl999l) . These authors 
compared the yields from core-collapse supernova models us- 
ing two different reaction rate libraries. For the 15 M Q stellar 
model studied, the two calculations gave 44 Ti yields that agreed 
to within 20%, in spite of the fact that many individual nuclear 
reaction rates differed by a factor of two or more between the 
two rate compilations. 

On the other hand. lNassar et alJ (l2005t) have recently mea- 
sured the 40 Ca(a', y) 44 Ti reaction cross section in the energy 
range for nucleosynthesis in supernovae. In that energy range, 
the authors find that the reaction rate is 5 - 10 times larger than 
the previously used theoretical rate calculated from a statisti- 
cal model ( Rauscher et al. 2000). This large difference between 
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the experimental rate and the theoretical rate may be due to the 
fact that the low level density in the 44 Ti compound nucleus 
limits the applicability of the statistical model for theoretical 
predictions for the rate of this reaction. In any event, the larger 
rate inc reases the yield of 4 4 Ti by a factor of ~ 2 in the stellar 
models tNassar et a j] 120051) explored. Such a large increase in 
the 40 Ca(o', y) 44 Ti reaction rate may allow normal core-collapse 
supernovae to account for the solar system's supply of 44 Ca; 
however, this result would worsen the discrepancy between the 
observed Galactic 44 Ti gamma-ray flux and our predictions. 

3.5.3. The supernova explosion model 

In core-collapse supernovae, 44 Ti production occurs by the 
alpha-rich freezeout near the mass cut. The location of the mass 
cut in the star will then certainly affect the 44 Ti yield. Also im- 
portant is the question of whether the simple-minded notion of 
a mass cut at a single radial shell in the star even makes sense in 
more realistic models that account for large-scale fluid motions 
behind the stalled supernova shock prior to the explosion and 
for stellar rotation. These more realistic models suggest that the 
material ejected from near the mass cut will in fact be a mix- 
ture of parcels that arose from both inside and outside the mass 
cut. We can certainly expect variations in the entropies of those 
parcels, w hich, in turn will h ave attendant variations in the 44 Ti 
yield (e.g, Pru et et all2005t) . 

For our purposes, the important issue to consider is how 
much variation can we expect in the 44 Ti yield from differ- 
ences in the mass cut and multi-dimensional effects. One- 
dimensional models suggest that typically half or more of the 
44 Ti produced during the explosion might fall back on the rem- 
nant (e.g, iNassar et alJl2005l) . Similar results are possible for 
the multi-dimensional models. It is therefore quite conceivable 
that yields of 44 Ti from supernovae of the same mass might 
vary by factors of > 2 simply due to variations in the location 
of the mass cut or multi-dimensional effects. Of course, if the 
supernova forms a black hole with mass greater than ~ 2 so- 
lar masses, it will swallow up its innermost material and, thus, 
most or all of the 44 Ti (and 56 Ni) it produced. Such supernovae 
would be dim in both visible and gamma radiation. Perhaps the 
Galactic 44 Ti map is indicating that supernovae over the last 
few hundred years have been predominantly of this type. 

3.6. Supernova Homogeneity 

Although often taken for granted, homogeneity among super- 
novae of a type remains an open issue: For thermonuclear 
supe rnovae, light c urves have been found to be fairly sim- 
ilar jBranchl Il998t) . Their successful empirical relative ad- 
justment through a light-curve-decline parameter apparently 
makes them "standard candles" over the full range of cos- 
mic evolution (This is the basis for the dete rmination of 
cosmic expansion history, see e.g. iDahlen et alJ §004)). The 
homogeneity of the r-process elemental abundance pattern 
in low-metallicity stars suggests that the r-process, which 
is commonly attributed to core collapse supernovae, also 
presents a fairly well-regulated nucleosynthesis environment 



(e.g. iThielemann et all l2002h . On the other hand, the 56 Ni 
masses ej ected in supernova e appear to scatter, within 30% 
for SNIa (Ben etti et al.l2 005). and f or core collapse events over 
a wide range from 0.01 to 1 M n dWoosley & Weaverlll995t 
IWooslev et alll99^lThielemann et alJr996l) . suggesting more 
variability in the core colla pse nucleosynthesis than in ther- 
monu clear explosions jNomoto et al.lll997t IThielemann et alJ 
2004). Still, for supernovae of type la alternative model 
types are also discussed that would produce quite distinctly 
different nucleosynthesis products than centr al carbon igni- 
tion in Chandrasekhar-mass white dwarfs jNomotol Il982t 
iLivne & ArnetJl995tlWooslev & Weaver! 1994 . 

In summary: supernovae of both types are quite homoge- 
neous in some of their characteristics, but anomalies suggest 
a deeper study of physical regulations and their observational 
impact. 

3.6.1 . 2-3D Effects in Core Collapse Supernovae 

Naga taki et af] dl997h in their explosive nucleosynthesis cal- 
culations of 2-D axisymmetric Type II supernova found that 
materials engulfed by energetic shock waves along polar di- 
rections undergo higher temperatures (or higher entropy per 
baryon JFrver & Heg erlboOO)) to produce a higher amount of 
^Ti th an in spherical explosions. Recently. iMaeda & Nomotol 
(2003) also studied hydrodynamics and explosive nucleosyn- 
thesis in bipolar supernova/hypernova explosions. Their bipo- 
lar models produce a large amount (>10~ 4 Mq) of 44 Ti and 
at the same time eject a relatively small amount (~0. 1-0.2 
Mq) of 56 Ni. These features of 2D supernova models in- 
spire IPrantzosI d2004l) to suggest that the "missing 44 Ti " 
problem (to be in concordance with SN1987A and Cas A 
observations and also to account for the 44 Ca solar abun- 
dance) could be solved by avoiding the overproduction of 
56 Ni. More systematic studies are therefore required under 
variety of progenitor masses, explosion energies, metallici- 
ties, and other physical variables. As this 2D effect could 
elude the discrepancy between the 44 Ti production in spher- 
ical supernova models and the amount observed from Cas 
A SNR in gamma-line fluxes and the amount deduced from 
44 Ca in solar abundance, the axisymmetric explosion also 
seems to be the natural consequences of rotation and mag netic 
field effects during pre-super nova phase ( Miiller & Hillebrandl 
Il98lt f^mada & Satolll994l) and the neutrino-driven convec- 
tion ( Bu rrows et al.l 1 19951) for the core-collapse supernovae. 
Furthermore, evidence to support the axisymmetric explosion 
have been inferred from the measurements of pulsar veloci- 
ties dBurrows & Haveslfl99d iHobbs et alJ l2004) and from the 
jet features observed in radio and X-ray im ages from several 
supernova remnants dGaensler et aljfl998l) . With various de- 
grees of Rayleigh-Taylor instabilities that may develop in core- 
collapse supernova (either in 2-D or 3-D simulations), we ex- 
pect there would be a wide distribution of 44 Ti production in 
supernova events as inferred indirectly by the pulsar velocity 
distribution. Thereby the task is to find 44 Ti distribution produc- 
tion with supernovae synthesizing higher than typical amount 
of 44 Ti produced in a spherical model (to explain the amount 
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of Ca in solar abundance and Cas A SNR), but also produces 
small amount of 44 Ti for the most-recent supernova so their 
gamma-line emissions are too weak to be detected. For this, 
we have to wait until the explosive 3-D nucleosynthesis can be 
performed within reasonable time. 

3.6.2. Rare 44 Ti rich events 

The apparent deficit of 44 Ti remnants from typical supernovae 
as the main source of 44 Ca is explained if typical supernovae 
are not the main source. Some rare type of event with a propor- 
tionately higher 44 Ti yield could be t he major source and no t 
leave detectable remnants today (e.g. JWooslev & Pintoll9 88). 
For example, these could be He-triggered sub-Chandrasekhar- 
mass thermonuclear supernovae JWooslev & Weaverlll994l) . If 
their recurrence time is several times the 44 Ti lifetime, the 
Poisson p robability of havi ng none detectable now can be large. 
IWooslev & Weave J (1 1994 find some models with up to one 
hundred times higher 44 Ti yields than the typical supernova 
values discussed above. These need only occur now every one 
to two thousand years to provide the necessary 44 Ca synthe- 
sis. Such a scenario, while potentially invisible to gamma-ray 
astronomy, would imply inhomogeneities among the relative 
abundances of 44 Ca and other isotopes (see section [D.3> . A 
galactic survey at 6.9 keV could test this scenario (Leising 
2001). Another a-rich freezeout nucleus, 59 Cu, decays to 59 Ni 
whose half-life is 75,000 years. The nearest remnants of these 
rare objects from the past 10 5 years could be easily detectable 
in the subsequent cobalt K a x-rays. 

Mainstream SiC grains, however, argue against rare pro- 
ducers of large amounts of 44 Ti being responsible for most of 
Galactic 44 Ca, as the 44 Ca/ 40 Ca ratio in mainstream SiC grains 
does not vary much from grain to grain; here sample size is 
much larger than for X grains, so a more representative sam- 
pling may be assumed, (see AppendixlDl. 

We perform a simple simulation to test the viability of the 
rare-event scenario. In this model, we take the model that pro- 
duces Fig. IB .41 (where the amount of 44 Ti in supernovae are 
as produced by the supernova models explained in Sec. 12. 51 
but modify the 10% of core-collapse supernovae to produce 
20x the 44 Ti of the supernova models. Therefore on average 
the amount of 44 Ti in this model is about a factor 2.7 of the 
model shown in Fig. IB .41 to account for the 44 Ca solar abun- 
dance. This model gives a better agreement with the observed 
44 Ti sky than the model used to produce Fig.[2]where all super- 
novae produce 3x the amount of 44 Ti of supernovae models. 
As shown in the caption of Fig.|2l the probability of this model 
detecting 0, 1, and 2 44 Ti sources are 0.012, 0.053, and 0.12, 
respectively. This results show that the model is more proba- 
ble than the model of Fig.[2]in explaining why COMPTEL and 
INTEGRAL only detect one 44 Ti source (less than expected). 
For comparison, the supernovae recurrence rates of the peak of 
the probability curve of the 1.157 MeV y-line fluxes that are 
consistent with the COMPTEL observed fluxes (see Appendix 
IB.2l and the dotted curve in Fig. IB.4t for the model of Fig. IB.4l 
this model, and the model of Fig. 13 are 36 yrs (Fig. IB.4L 40 
yrs, and 58 yrs, respectively. This result is encouraging that a 



set of supernovae models with a rare type supernova that pro- 
duces most of the solar 44 Ca gives a closer recurrence rate to 
the one implied by the historical supernova record (~17 yrs) 
than the model with 3x amount of 44 Ti of the calculated super- 
nova models. Still we have no good explanation for why some 
small number of supernovae are so rich in 44 Ti . 

4. Conclusions 

The observed distribution of pointlike sources of 44 Ti in our 
Galaxy is inconsistent (at a probability near 10~ 3 ) with the 
combined current understanding of 44 Ti production, specifi- 
cally 

- core collapse supernovae eject 10~ 4 M of 44 Ti , within a 
factor of two, and 

- the galactic core collapse SN rate is about 3 per century, 

assuming the last few centuries are representative of the steady- 
state 44 Ti production. The inconsistency is exacerbated if we 
further demand, as currently understood to be so, that the stan- 
dard solar abundance of 44 Ca originates from 44 Ti -producing 
supernovae. The larger discrepancy is because the product of 
the above two quantities falls short by a factor of three of 
the requisite current 44 Ca production rate. The disagreement 
persists for any combination of rate and yield whose prod- 
uct is the required value, unless the yield of 44 Ti is very 
much higher and the rates very much lower, i.e., 44 Ca is not 
made primarily by typical supernovae events, but by very rare 
ones. These might include He-triggered detonations of sub- 
Chandrasekhar SNe la, or rare variants of core collapses, per- 
haps those most departing from spherical symmetry. Evidence 
of He-cap SNIa as source of the 44 Ca abundance might be iden- 
tifiable in 44 Ca/ 40 Ca ratios greater than solar in some main- 
stream SiC grains. 

A future survey with gamma-ray line sensitivity of 
10~ 6 cmT 2 s~ l would be expected to detect > 10 sources (Fig. 2), 
and so could rule out 44 Ca production by frequent supernovae 
at confidence 5 10 . Regardless of specific assumptions, the 
Cas A supernova remnant as the brightest 44 Ti remnant in the 
galaxy is a priori very unlikely. That the brightest SNR should 
be found in the outer galaxy or that it is over 300 years old are 
each improbable. It suggests that yields higher than suggested 
by many current calculations are possible, which, of course, 
makes the lack of other detectable remnants even more puz- 
zling. 
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Appendix A: Simulations of 44 Ti Skies 



A.2. Supernova Types and Spatial Distribution 



A.1. Monte Carlo Simulation Details 

In generating supernova events of our Mont e Carlo simulations , 
we adopt the procedure shown by | Higdon & Fowled (^987); 
Mah onev et alJdl992l) : lHartmann et alTdl993R . A random num- 
ber between and 1000 is uniformly generated to represent a 
supernova age between zero and 1000 y 1 1 .4xt7-,44) where 
ttm = (87.7 + 1.7) yrs dNorman et alJl997HGorres etall 19981 
lAhmad et all 1998) is the mean life of ^Ti. This range is large 
enough that a supernova age older than 1000 y does not con- 
tribute to the Galactic 44 Ti flux. 

Another random number is generated to choose the type of 
the supernova event. Then, several random numbers are gen- 
erated to give us the location of the supernova according to 
its spatial distribution. A detail procedure in generating the lo- 
cations of supernova from a disk and spheroid populatio ns of 
Type la events can be f ound in the paper of Higdon & Fowler 
( 1987k iMahonev et all dl992l) . The distance of the supernova 
location to the Sun can be calculated easily. 

A random number is generated from a Gaussian distribu- 
tion to give us the peak bolometric magnitude of the supernova 
following the distribution given in section |A~21 The apparent 
magnitude then can be calculated knowing the location, the 
distance, extinction magnitude (from Hakkilla et al.'s empirical 
model), and the peak bolometric magnitude of the supernova. 
We find in model A (Section IA.2> the fraction of SNe in the 
last millennium brighter than apparent magnitude is ~11% 
(see Fig. |B3}. 

The amount of 44 Ti of the supernova is given by a uniform 
distribution of 44 Ti according to its type as described in section 
12.51 however note that the analysis in this appendix does not 
use the multiplying factor of 3 that is used in \2.5\ F y , the 44 Ti 
gamma-line flux at the Sun location, then can be determined 
knowing the location, the distance, the age, and the amount of 
44 Ti of the supernova: 

F y = 8.21 x 10~ 3 M 4 exp(-t/87.7yr)/d 2 kpc y cmT 1 s~ x 
where M\ is the 44 Ti of the supernova in units of 10~ 4 Mq and 
dkpc is the distance (in kpc) of the supernova from the Sun. We 
find that in Model A the fraction of supernovae in the simula- 
tions with 44 Ti y-line fluxes larger than lxl0~ 5 y cmT 2 s~ l and 
3xl0~ 5 y cmT 2 s are ~11% and ~5%, respectively (Fig. [3}- 
Its longitude distribution along the Galactic plane in some rel- 
evant 44 Ti gamma-line flux bandwidths can be seen in Fig.0] 
while for models B and C (Section lA.2> are shown in Fig. IA.ll 

We generate 1 million supernovae according to the above 
prescription and construct galaxies with numbers of super- 
novae between 1 and 300. For each galaxy with a certain n 
number of supernovae, the 44 Ti gamma-line flux distribution 
F(fy,n) and the SN apparent magnitude distributed M(m,n) 
can be extracted. 

For a galaxy with SN recurrence time of T rec or average su- 
pernova rate, p = 1000/T rfc (we only simulate supernova events 
in the last millennium as described above), the 44 Ti gamma- 
line flux distribution, F(f y ,p) can be obtained from F(f y ,p) = 



Supernova events as the source of Galactic 44 Ti can either be 
Type la, Ib/c, or II SNe. For this study, we choose a generally 
accepted (though uncer tain) value of the type ratio, la : Ib/c : 
0. 1 : 0.15 : 0.75 dPawson & Johnsorj|l994l lHatano etal 



II 



ll997t) . lYoshiietalJ ( 1996) and Tsuiimoto et al 



( 119951) using 

their chemical evolution model find that the ratio of the to- 
tal number of Type la to Type II SNe of 0.12 gives a good 
agreement with the observed solar abundance. This is consis- 
tent with the observed la frequen cy which is as low as 10% o f 
the total SNe occurrence Jvan den Bergh & Tammanr]ll99ll) . 
Monte Carlo representations of Type la SNe are generated us- 
ing a nova distrib ution template that trace s the blue light dis- 
tribution in M31 llHigdon & Fowlerlll987l) . These populations 
form an axisymmetric disc and a spherically symmetric bulge. 
The spheroid density distribution follows: 

„ ,.. -10.093[(/-)" 4 -l] 

- 1.25 (^r 6/8 e R © for R < 0.03% 

-10 093[(-£-l" 4 -l] 
_ fJLV 7 / 8 P U «Q ; J r1 _ 0.08669 

*0 



[l-T^]for»>0.03fi o 



and the disc density distribution follows: 

- n(p,z) oc e - w<r * e~ {p - R Q )/Ph 

- cr z = 0.325 kpc (Mihal as&Binnevll98ll) 

- pi, = 3.5 kpc (de Vaucouleurs & Pence! 1978 



The fraction of SNIa o ccurring in the spheroid is taken to be 
-1/7 of the total SNIa feahcall & Soneirall980l) . 

Supernovae of Types lb & II are associated with massive 
stars whose birth places are exponentially distributed in height 
above the plane with a scale length of ~100 pc. Because the 
rate of core collapses is larger than the rate of thermonuclear 
supernovae, we study several Type II distributions to ensure 
that our results do not depend significantly on this choice. We 
consider three cases: 

1 . Model A: Exponential disk with no supernova within 3 kpc 
of the Galac t ic cen ter, with radial scale length of 5 kpc 
lHatano etal Jll997t: 

n(p, z) oc g-W/°v e -p/pi 
cr z = 0.100 kpc 
Pd = 5.0 kpc 

2. Model B: Exponential disk with radial scale length of 3.5 
kpc that pr oduces an acceptab le fit to the COMPTEL's 26 Al 
y-line map Diehl et al] ll 19951) . 

3. Model C: Gaussian-ring disk at radial distance of 
3.7 kpc and radial di stance scale length of 1.27 kpc 
llTavlor&Cordesll993l) . 

The distribution of supernovae in model A in directional is 
shown in Fig. lA.2l and in distance from the Sun is shown in Fig. 
lA3l 



Appendix B: Comparison of Supernova Rates from 
44 Ti Gamma-Rays and the Historical Record 

In order to produce a consistent picture of the recent 
Galactic supernova rate with the best known 44 Ti supernova 
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Fig. A.3. The distribution of supernovae distance in model A. 



yields, supernovae types and spatial distribution, and Galactic 
extinction model we utilize our Monte Carlo simulation 
method to compare two observables of the same phenomenon 



faigdon & FowleJl987tlMahonev et alJl992tlHartmann et all 
Il993t iThe et alJll999l) . We use the COMPTEL gamma-ray 
maps and the six historical Galactic supernova record of the 
last millennium to constrain the range of Galactic supernova 
rates. 

Similar to the above determination of 44 Ti flux distribution 
from a simulated set of supernova events, other distribution 
such as the apparent magnitude distribution of supernovae for 
a galaxy with an average supernova rate p can be constructed 
(Section IB.3> . From this distribution we can determine the 
probability of optical detection for a certain detection-limiting 
apparent magnitude, such as magnitude 0, for example. In this 
way, we can use the same underlying assumptions about su- 
pernova characteristics, and compare two different observables, 
i.e. 44 Ti gamma-ray versus optical detections of the supernovae 
at the adopted rates. Varying the rates then so that they conform 
to the observational constraints, we obtain a handle on system- 
atic differences in observed supernova rates through these two 
observational windows. 

For the 44 Ti gamma-ray observations, we utilize the 
COMPTEL survey discussed above (Sections l2. ll and lTTt . For 
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Table B.l. Recent Galactic Supernova Record 



Name 


Year 


Distance 


1 


D 


Type 






(kpc) 








Lupus (SN 1006) 


1006 


2.2 


327.57 


14.57 


la 


Crab 


1054 


2.0 


184.55 


-5.79 


11 


3C58 (SN 1181) 


1181 


2.6 


130.73 


3.07 


11 


Tycho 


1572 


2.4 


120.09 


1.42 


la" 


Kepler 


1604 


4.2 


4.53 


6.82 


Ib/II b 


Cas A 


1680 


2.92 


111.73 


-2.13 


lb 



*; Studies of X-ray emission from Tycho 's shoc ked eject a show the 
SNR was created by a Type la SN. \Hwane et alil99& \B adene s et all 

\200A\200A) 

b : The uncertainty of Kepler's SN Type was reviewed bv \BlaiA tooi) . 



(with which they are consistent.) For each Monte Carlo Galaxy, 
aj 2 value is calculated by comparing the two strongest fluxes 
with the fluxes measured from Cas A and GRO J0852-4642 and 
the other fluxes are compared with null fluxes. The probability 
of the model to be consistent with the COMPTEL fluxes is de- 
termined from the y 1 and using the number of degrees of free- 
dom as the number of independent 10°xl0° image elements 
(fields of view) of the supernova positions in the model. This 
size of independent fields of view is obtained by calibrating the 
size to produce the results of the analysis of COMPTEL's map 
of the inner Galaxy. Averaging the probabilities of 10 6 Monte 
Carlo galaxies, we obtain the average probability shown in Fig. 
IB.4l and lB31 as dotted lines. 

B.3. Historical Record Analysis 



the optical observations, we utilize the historical record sum- 
marized in Table lB.il 

For comparison of the COMPTEL gamma-ray map with 
the results of Monte Carlo simulations, due to our limited com- 
putational power, for the purpose to reduce systematic effects 
from regions of low exposure, and due to our limited ability to 
analyze the map, we perform two analysis. In the flux datas- 
pace (Section IB.2I we use the flux information exclusively, 
without using the observed location of the gamma-ray source. 
In the map dataspace analysis (Section lB.il we only use the 
COMPTEL map of inner galaxy (|Z| <90°, \b\ <30°) which does 
not include the locations of the detected 44 Ti y-line from Cas A 
and GRO J0852-4642. In order to see how consistent the model 
with the Galactic supernova record, we perform the historical 
record analysis below. 

B.1. Map Dataspace Analysis 

In map (imaging) analysis, a simulated COMPTEL data set 
is produced by convolving the directions and the 44 Ti y-line 
fluxes of the Monte-Carlo-generated supernovae through the 
COMPTEL detector response. The probability of the consis- 
tency of the simulated data with the measured counts is car- 
ried out by calculating the likelihood of the simulated data plus 
the best-fit background model, and comparing t his to the l ikeli- 
hood of the background model only Jde Boer et al.ll9 92). The 
change in the likelihood gives the relative probability that the 
specific realization of the model is consistent with the mea- 
sured COMPTEL data. This process is repeated for at least 10 4 
galaxies to obtain the average relative probability of the model 
at a particular supernova rate as dashed lines in Fig. IB .4! and 

1531 



la : lb 



0.10 : 0.15 : 0.75 



B.2. Flux Dataspace Analysis 

In flux dataspace analysis, in order to include the detected 
sources, the locations of the 44 Ti supernovae are not used, flux 
information is used exclusively. The flux dataspace consists of 
the t wo 44 Ti detec t ed flux es from Cas A SNR and GRO J0852- 
4642 Iyudin et al. ( 1998) with other fluxes assumed to be zero 




Fig. B.l. The distribution of the peak absolute magnitudes of 
supernova in the B-band of model A. 



The historical supernova record covering the last millen- 
nium shows a total of only six Galactic SNe during that era 
(Table lB~TT i. Of course, this small number is due to significant 
losses from extinction and incomplete monitoring of the sky, 
especially during the early centuries. In fitting this data, we 
count the fraction of galaxies that have six SNe brighter than 
magnitude 0. In this approach, we assume that historical SNe 
were detected if they were brighter than magnitude 0. 

The peak absolute magnitudes of supernovae in the B-band 
are approximated by Gaussian distributions. For the mean val- 
ues and the one standard deviation we adopt for type la Mb = 
- 19.4, cr=0.2 (lBranchlfT 998). for type lb M B = -18.2, cr=0.3 

1 Dawson & Johnson! 1994TTand fo r type II M B = -17.2, cr=1.2 
van den Bergh & Tammannll99ll) . Their distribution is shown 
in Fig. IB. II Observed magnitudes of simulated SNe are ob- 
tained by convolving absolute magnitudes through a Galactic 
extinction model of llHakkila et alJll997l) . In this extinction 
model, the total visual extinction from the Sun's location to 
the Galactic Center is 1 1.64 magnitude and the total visual ex- 
tinction perpendicular to the Galactic plane at Sun's location 
is ~0.1 magnitude (Fig. IB .2b . The cumulative distribution of 
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supernovae having blue peak apparent magnitude brighter than 
oig in model A is shown in Fig. IB. 31 



Ib/ll SNe: Exponential Disk for R > 3 kpc 
la : lb, : II = 0.10 : 0.15 : 0.75 , 
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Fig. B.2. 77ie Galactic visual extinction magnitude as a func- 
tion of distance fro m the Sun to various direction according to 
empirical model oflHak kila et al\ 1997). 
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Fig. B.4. Normalized relative probability of Galactic SNe rate 
calculated from three different analysis vs. SNe recurrence time 
for model A. The solid line is the probability of the model in 
which the number of SNe with my < is 6 in a 1000 year 
duration. The dashed-line is the probability inferred by maxi- 
mum likelihood analysis of COMPTEL 's 1.157 MeV y -line in- 
ner galactic (I <90, b <30) image map. The dotted-line is the 
relative probability from chi-square test that the 1.157 MeV 
y-line fluxes in the model are consistent with the COMPTEL 
observed fluxes. The labels of the curves show the optimal, 
Icr ranges, or the lo~ lower limits of the SNe recurrence time. 
The boxes, circles, and diamonds on the curves show the lcr, 
90%, and 95% relative probabilities. The shaded area shows 
the range of the Galactic SNe recurrence time that is consis- 
tent with better than lcr probability with each of the historical 
record and the COMPTEL flux dataspace analysis. 



Fig. B.3. The apparent magnitude of superno vae in m odel 
A acc ording to empirical extinction model of lHakkila et al\ 
\1997\> . 



The average fraction of the model that is consistent with six 
events brighter than magnitude is shown in Fig. IB .41 and lB.5l 
as solid lines. 



B.4. The Supernova Rates: Optical versus 
Gamma-Ray Constraints 

Our Galactic supernova record analysis leads to a most proba- 
ble Galactic supernova recurrence time of ~ 17, ~ 16, and ~13 y 
based on models A, B, and C, respectively. This implies a larger 
rate than given by previo us investigations based on the histori- 
cal record. For example, Dawson & Johnson! dl994 estimated 
~3 SNe per century as also obtained by Tammann et al. ( 1994). 
The rate is smaller than ours because Dawson & Johnson 



(1994) considered 7 observed SNe within the last 2000 yrs and 
assumed that the historical record is 80% compl ete. A hi gher 
SNe ra te of 5 SNe per century was obtained by lHatano et alJ 
who include a population of "ultradim" SNe in addition 
to 4 observed Galactic SNe having V < and 80% complete- 
ness within the last millenium. 

The map and flux methods of estimating supernova rates 
based on y-ray data appear to be significantly different, but they 
are statistically consistent with each other for a wide range of 
supernova rates. The imaging analysis gives a smaller rate than 
the flux analysis, because the data used in the map analysis 
do not include the Cas A and GRO J0852-4642 y-line detec- 
tions. However, the probability estimate from the flux analysis 
is a somewhat coarser estimate because the size of independent 
FOVs used in determining the probability is not known exactly. 
Ignoring the spatial information begs the question as to why the 
two brightest SN are in the outer Galaxy rather than in the inner 
Galaxy where we expected them. 

In the flux dataspace analysis, where we ignore the ex- 
pected spatial distribution of 44 Ti remnants and consider only 
the measured flux distribution, we find a most probable super- 
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nova rate that is more compatible with standard values (such as 
the rate inferred by the historical record), even for standard 44 Ti 
yields. However, the COMPTEL map data indicate a lower SN 
rate (~1 SN/36 yr) than that suggested by the historical record 
(~1 SN/17 yr). However model A which does not have a su- 
pernova event in the inner 3 kpc radius has better agreement 
between the map analysis and the historical record analysis. 
This suggests that a model with less concentrated supernova 
near the Galactic center than the model used here could give 
a better agreement than we obtained here. Also, an extinction 
model with smaller extinction toward the Galactic center than 
Hakkilla et al.'s model could give a better agreement than what 
we present here. 

Based on chemical evolution studies, Tim mes et alJ |l996) 
estimated that only ~l/3 of the solar 44 Ca abundance is ac- 
counted for. Models with a SN rate of ~3 SNe per century and 
standard 44 Ti yields fail to produce the solar 44 Ca abundance. 
This rate, when confronted with the gamma-ray data (dashed 
line in Fig. IB.4l and lB.5l l is too large: the COMPTE L gamma- 
ray da ta worsen an already serious problem. iTimmes et alJ 
(1996) suggest 3 possibilities: 

1. Increase the 44 Ti yields by a factor of ~3. 

2. Increase the supernova rate by a factor of ~3. 

3. There is another source of 44 Ca in the Galaxy. 

Our analysis in Fig. IB.4l and lB.5l shows that the first and sec- 
ond option are not compatible with COMPTEL's 44 Ti y-line 
map, which would be brighter by a factor of 3 or exhibit a 
larger number of 44 Ti hot spots than actually observed. We are 
thus left with the third option, to seriously entertain the idea 
that there exists some rare type of supernova (i.e., detonation 
of helium white dwarf), not realized in recent centuries, that 
produces very large amounts of 44 Ti. 

Appendix C: Solar 44 Ca constraints. 

We now wish to estimate the effects of choices of chemical evo- 
lution model parameters o n the 4 4 Ca constraint. We do thi s with 
analytic models faavtonlll985l) .lL eising & Shard ill 994 used 
that chemical evolution model to estimate the current 44 Ca pro- 
duction rate by requiring the yield that gave exactly the solar 
abundance 4.55 Gyr ago. With model parameters k (infall con- 
tribution) and A (time delay) defining the shape of the galactic 
infall history, the current production is 

pCCd) = X ( 44 C«) *±i M g (T G ) [Zk±^_ ( Zk±Vri, (C. 

where T G is the age of the Galaxy, T Q =T G -4.55 Gyr, and 
M?(T G ) is mass of interstellar gas participating in star for- 
mation now. In this formulation, the 44 Ca production depends 
mainly on the infall parameter k= [f(t)/M G (t)] (t + A) and the 
gas mass, and weakly on A and the age of the Galaxy. The 
closed-box model, k=0, is ruled out by observati ons, and k — 2- 
4 is f avored from a number of considerations dClavton et alJ 
Il993l) . For the total 44 Ca production rate, therefore, a low 
extreme value of p( Ca)=l.l 1CT 6 M yr -1 is obtained for 
k=l, A=1.0 Gyr, M c (r c )=4 1O 9 M , and T c =13 Gyr. An up- 
per limit is estimated as p( 44 Ca)=1.2 1O~ 5 M yr -1 for k=4, 



A=0.1 Gyr, M G (T G )=1 10 I0 M o , and T c =10 Gyr. Our favored 
value is p( 44 Ca)=5.5 1(T 6 M yr" 1 , derived for k=2, A=0. 1 Gyr, 
M G (r G ) =l lO lo M ,andTp=12 Gyr. As in previous such stud- 
ies (e.g. Leising & Sharel fl994) we assume that this equals 
the 44 Ti production rate. (However, we note the statement by 
ITimmes et alJ i 19961) that one-half of ^Ca is made directly as 
44 Ca rather than as 44 Ti.) In such models, the solar 44 Ca abun- 
dance thus could be provided by, for example, 2.8 SNe per cen- 
tury at the present epoch ejecting on average 2 1O~ 4 M of 44 Ti. 

Appendix D: Presolar grains 

Do the presolar SiC grains dClavton & NitMerl l2004h 
that condensed within the expanding s upernova interior 
jAmari & Zinnerll997tElavton et all 1997 a) shed light on the 
44 Ti origin? Called "X grains", each represents a sample of 
selected portions of a supernova interior, but "what supernova" 
and " what portions of it" are unidentified questions for each X 
grain. iDeneault et alJ ll2003l) present a physical condensation 
argument to identify the region within the supernovae; but 
their argument awaits digestion by the condensation-chemist 
community. The X grains certainly contain evidence of having 
condensed with abunda nt five 44 Ti in t he form of very l arge 
44 Ca/ 40 Ca ratios JNittleretal J 11994 IClavton et alJ Il997at 
IHoppe et alJ l200(il) that can be explained only through the 
action of live 44 Ti. This proves their supernova origin, and 
was a predicted signa ture of supernova origin for presolar 
grains dClavtonl [1,975) from the beginning. If each grain 
contains sufficient Ca and Ti the SIMS analysis can yield the 
^Ti/ ^Ti abundance ratio at the time of conde nsation ( Fig. 
2 of IClavton et alJ Jl997al) and Figs. 8 and 9 of IHoppe etall 
feOOd) ). Many X grains having 44 Ti/ 48 Ti = 0.1-0.6 are found. 
The production ratio required by the assumption that all 44 Ca 
within solar abundances is the product of decay of 44 Ti is 
P( 44 Ti)/P( 48 Ti) = 0.72; therefore, supernova X grains are found 
with ratios up to the required bulk production consistent with 
that assumption. Despite this it must be clearly understood 
that the measured production ratio P( 44 Ti)/P( 48 Ti) within 
material in a specific grain bears no simple relationship to 
the bulk production ratio in the supernova within which that 
grain condensed. This illustrates the problem of not knowing 
precisely what supernova material the grain's isotopes reflect. 
However, this proble m seems capable of eventual solution 
JPeneault et all2003l) . 

b.1. SNIIX grains 

The SiC X grains probably originated in core- collapse S NII. 
Withi n them the pure alpha-rich freezeout dWooslev et alJ 
Il973l Section VII), which produces the bulk of ^Ti ejected 
from SNII, has production ratio near P( 44 Ti) /P( 48 Ti) = 1-2 (see 
Fig. 23 of IWooslev etai]ll973l iThe et alJll998l Section 3); 
so no X grain is pure alpha-rich freezeout material, although 
some of them must contain a significant fraction of their Ti 
from the alpha-rich freezeout. On the other hand, the produc- 
tion ratio during normal silicon burni ng is very mu ch smaller, 
near P( 44 Ti)/P( 48 Ti) = 0.01 dWooslev et alJll973l Fig. 19). 
Therefore, X grains having 44 Ti/ 48 Ti near 0.01 or less may 
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contain no alpha-rich freezeout material at all, but instead con- 
dense from Si-burning ejecta rich in 28 Si. 

In t heir study of 99 SiC X grains from the Murchison me- 
teorite, iHoppe etalJ J2000) found that 25 contained enough Ti 
and Ca for isotopic analysis with the CAMECA IMS3f ion mi- 
croprobe in Bern. Now that the new nano-SIMS ion microprobe 
is functional, we may hope for even more complete future sur- 
veys with its higher sensitivity. Of the 25 having enough Ti and 
Ca, 5 revealed large and easily resolved 44 Ca excesses, corre- 
sponding to P( 44 Ti)/P( 48 Ti) > 0.01. Although the chemistry of 
conden sation of SiC in supernovae has not been solved, initial 
studies ( Dene ault et alJ l2003) of the location of C and Si sug- 
gest that no more than 20% of SiC X grains should be expected 
to contain P( 44 Ti)/P( 48 Ti) ratios in excess of those available 
during normal O and Si burning. That is, a majority of the X 
grains should contain no alpha-rich matter even if each super- 
nova ejects such matter. If this be taken as so, the data suggest 
that the parent supernova population responsible for the 100s 
of X grains that have been studied did eject alpha-rich freeze- 
out material in addition to its Si-burning matter. If ejection of 
alpha-rich freezeout matter were instead rare, we expect that 
the number of SiC grains having large P( 44 Ti)/P( 48 Ti) ratios 
would be significantly less than 20%. Without more detailed 
understanding of the origin of the X grains it seems plausible 
to take their evidence to suggest that all SNII eject alpha-rich- 
freezeout matter. 

If the SiC X grains do originate in SNII, their frequent 44 Ca 
excesses demonstrate that either the X grains all originated in 
a single nearly supernova that did eject 44 Ti or originated in 
many presolar SNII that mostly ejected 44 Ti (although not nec- 
essarily always). The properties of the "mainstream SiC grains" 
enable us to argue against the possibility of a single supernova 
as origin of the X grains. The much more abundant mainstream 
grains, which are thought to have originated in presolar carbon 
stars, have been argued on the basis of their Si isotopes to have 
been the result of a great many C stars. Their Si isotopic com- 



positions represent approximately the initial Si isotopic com- 
positions from which the intermediate-mass stars formed, and 
their correlation between excess 29 Si and 30 Si is explained by 
galactic chemical evolution effects that have produced a wide 
range of Si isotopic compositions in a larg e number of C stars 
( Timm es & Clavtonlll996t IClavtonlEo03l Because the main- 
stream grains comprise about 10~ 4 of all interstellar Si, much 
of which is quite old, and because the condensed SiC from the 
C-star phase is but a small fraction of all Si ejected from stars, 
it seems that the lifetime of SiC grains in the ISM is not short. 
If this be so, the lifetime of X SiC grains is also not short, sug- 
gesting a substantial number of contributing supernovae to the 
matter gathered into any cloud in the ISM. 

D.2. SNIa X grains 

For the question of 44 Ti gamma-ray hotspots the issue arises 
whether the SiC X grains may instead have arisen from ex- 
ploding SNIa. After all, the high expectation for seeing ^Ti 
supernova remnants shown in Figure 1 came from the assump- 
tion that Type II are the sources of 44 Ti nucleosynthesis and 
therefore also of the solar 44 Ca abundance. The absence of such 
sources could be explained if rare Type I events produced much 
larger yields than are expect ed from Type II. Do the supernova 
grains offer any guidance ? IClavton et al.l dl997bl) showed that 
isotopically good fits to X grains might originate within He 
caps on exploding CO white dwarfs. Huge production ratios, 
up to P( 44 Ti)/P( 48 Ti)=100, occur in those He shell zones hav- 
ing peak Ty >1 (see their Fig. 6). These are too great for exist- 
ing observations of X grains, suggestin g that SNIa are not their 
sources. However. Idavton et alJ dl99 7b) showed that good fits 
to all isotopes require that the He caps undergo considerable 
post-explosive mixing prior to condensation of the X grains. 
Their Table 2 shows production ratios as small as 0.01 in an 
average of zones 1-8 (the coolest eight zones) and near 10 in an 
average over all sixteen zones. (The reader must note that much 
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Ti production is listed by them under Cr, its radioactive pro- 
genitor.) The isotopic possibility must therefore be addressed 
that the 44 Ti-rich X grains are from SNIa He caps rather than 
from SNII. The discovery of even a single X-type SiC grain 
containing an initial ratio 44 Ti/ 48 Ti=5 or greater would demon- 
strate that the He-cap SNIa do exist and that they can condense 
SiC. This would amount to an existence proof for these rare 
44 Ti producers. But no such grain has yet been detected. 

IClavtonet"ai1lll997bl) pointed out that if the SNIa are near 
the Chandrasekhar mass, the He cap can be no more than 
0.01 M so that less than 10~ 5 M o of 44 Ti is ejected. Such 
events would be both rare and dim in 44 Ti lines. However, 
a much larger He cap is involved if the white dwarf is sub- 
Chandrasekhar (larger) and the detonation begins in a massive 
He cap. Large amounts of 44 Ti are ejected from such models 
JWooslev & Weaverll994 . so that the rates of such events are 
limited by the requirement that they not overproduce the galac- 
tic 44 Ca abundance. If this occurs, the rarity of these events in 
time could account for their absence in the COMPTEL data 
(see also the Discussion). The relevance of the X grains derives 
from the physical unlikelihood of condensation of X grains in 
such He caps. Their expansion is very fast, and the radiation 
environment intense, so that even the formation of molecules 
seems unlikely until the density has b ecome too low for the 
growth of a 1 -micrometer grain of SiC. IClavton et alJ Jl997b) 
note that even the possibility of grain condensation would seem 
to require 3-D explosive modeling in order that some of the 
He cap can remain at low velocity ("slow He"). This entire 
problem will require more study before plausibility of SiC X 
grains from SNIa can be admitted. We therefore conclude that 
the 44 Ti-bearing X grains have arisen in SNII. 

D.3. ISM Inhomogeneity of u Ca from He-cap SNIa 

We turn now to non-supernova Stardust whose 44 Ca ratios re- 
flect the initial compositions of their donor stars. The super- 
nova X grains seem to support the assumptions that led to the 
conflict between the COMPTEL map (Fig.[Q and the expected 
hotspots (Fig. |2j; namely, that X grains condensed within SNII, 
that the number of contributing SNII was large, and that most 
ejected alpha-rich freezeout matter. However, X grains do not 
demonstrate that the bulk ejecta of SNII contain sufficient mass 
of 44 Ti to account for the natural 44 Ca abunda nce. To test for 
the pr esence of rare large 44 Ti producers dWooslev & Weaver! 
1994) responsible for making good the shortfall from SNII we 
turn to the mainstream SiC grains with the following original 
argument. Rare SNIa responsible for roughly 2/3 of galactic 
44 Ca plausibly result in Galactic inhomogeneities in the inter- 
stellar 44 Ca/ 40 Ca ratio that should be larger than those seen 
by astronomers for elemental-ratio variations attributed to in- 
homogeneous incorporation of SNII ejecta. It must be slower 
to homogenize the ISM if SNIa occurring every 3000 yr con- 
tribute 2/3 of 44 Ca abundance than for homogenizing ISM from 
SNII occurring every 30 yr and making but 1/3 of the 44 Ca 
abundance. ISM regions temporarily enriched in the rare SNIa 
ejecta may make stars having larger 44 Ca/ 40 Ca initial ratios. 
The mainstream SiC grains reflect the interstellar composition 



from which the carbon stars formed. That is, factor of two or 
more variations in 44 Ca/ 40 Ca ratio should be present in the ini- 
tial compositions of stars. 

It is this expectation that th e mainstream grains may 
speak against. Ho ppe et alJ (2000) measured Ca isotopes in 
28 mainstream SiC grains (see their Fig. 8) and found their 
44 Ca/ 4l) Ca ratios to be indistinguishable from solar despite 
large Si isotope variations in the same stars demonstrating 
that the grains come from different AGB C stars having dis- 
tinct chemical evolu tion histories for their initial compositions 
jTimmes & Clavtonl 19961) . Each mainstream grain was consis- 
tent with solar 44 Ca/ 40 Ca ratio and their average was enriched 
3.8% in ^Ca, as expected for modest s-process enrichments of 
AGB atmospheres by the third dredge ups. Because of the small 
numbers of Ca atoms in mainstream grains, however, individ- 
ual grain ratios were uncertain by at least 20%. Even so, 20% 
is a small inhomogeneity for rare events producing 2/3 of the 
bulk 44 Ca. To improve this data base we call for high precision 
measurements of the 44 Ca/ 40 Ca ratio in a sample of mainstream 
grains with the new nanoSIMS ion microprobes. Such a survey 
could reveal more precisely any variations in initial 44 Ca/ 40 Ca 
ratios in intermediate-mass stars expected if the 44 Ti synthesis 
has instead been the result of a few rare events of very large 44 Ti 
yield. The present data speak against that. The new arguments 
that we have presented here can all be improved dramatically 
in the next few years. Although their message today is not sta- 
tistically certain, their sense is to support the original conflict 
that we have displayed in Fig. |2 Because the 44 Ca/ 40 Ca ratio 
in mainstream grains from many AGB stars is very near to the 
solar ratio, we point out that it is no longer possible to entertain 
the idea that the sun may itself contain an anomalous 44 Ca/ 40 Ca 
ratio. One can not blame the inability of standard models of the 
chemical evolution of the galaxy to produce enough 44 Ca by 
speculating that the sun is abnormally rich in its 44 Ca content. 

It should be noted that a thorough study of Stardust isotopic 
inhomogeneity within the context of inho mogeneous gal actic 
chemical evolution has been presented by Nittler (2005). He 
shows in particular that correlations between Si and Ti isotopic 
compositions found in Stardust can not be accounted for by in- 
homogeneous GCE. He does not, however, discuss the varia- 
tions in 44 Ca/ 40 Ca that are generated by his model owing to the 
He-cap SNIa. On the other hand, his calculations may not be 
ideal for this problem because his model admixes into a ho- 
mogeneous ISM the same fraction for each supernova's ejecta 
(a= 1.7xl0~ 5 ) after Monte Carlo sampling to obtain N=70 su- 
pernovae. Nonetheless, that paper shows the potential power of 
Stardust to delimit GCE inhomogeneities in general. This ap- 
proach should be reconsidered carefully when a good data set 
for initial 44 Ca/ 40 Ca in Stardust is available. 

Mainstream SiC grains are probably not the best grains 
with which to measure initial stellar 44 Ca/ 40 Ca abundance ra- 
tios. The mainstream SiC grains contain very little Ca because 
Ca is much less favored chemically than is Ti within thermally 
cond ensing SiC grains . The hibonite Stardust grains described 
bv lNittler et alJ J2005h . on the other hand, contain much more 
Ca because Ca is an essential ingredient of the hibonite crys- 
tal structure. Therefore, accurate measurements of 44 Ca/ 40 C in 
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hibonite grains could reveal better information relevant to ISM 
inhomogeneity of the 44 Ca/ 40 Ca ratio. 



